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ABSTRACT 

The spectral shapes of the contributions of different classes of unresolved gamma-ray emitters can 
provide insight into their relative contributions to the extragalactic gamma-ray background (EGB) and 
the natures of their spectra at GeV energies. We calculate the spectral shapes of the contributions to 
the EGB arising from BL Lacertae type objects (BL Lacs) and flat-spectrum radio quasars (FSRQs) 
assuming blazar spectra can be described as broken power laws. We fit the resulting total blazar 
spectral shape to the Fermi Large Area Telescope measurements of the EGB, finding that the best-fit 
shape reproduces well the shape of the Fermi EGB for various break scenarios. We conclude that 
a scenario in which the contribution of blazars is dominant cannot be excluded on spectral grounds 
alone, even if spectral breaks are shown to be common among Fermi blazars. We also find that while 
the observation of a featureless (within uncertainties) power-law EGB spectrum by Fermi does not 
necessarily imply a single class of contributing unresolved sources with featureless individual spectra, 
such an observation and the collective spectra of the separate contributing populations determine 
the ratios of their contributions. As such, a comparison with studies including blazar gamma-ray 
luminosity functions could have profound implications for the blazar contribution to the EGB, blazar 
evolution, and blazar gamma-ray spectra and emission. 

Subject headings: galaxies: active - gamma rays: diffuse background - gamma rays: galaxies 



1. INTRODUCTION 

The gamma-ray sky as currently observed by the Large 
Area Telescope (LAT) on board the Fermi Gamma- 
Ray Space Telescope consists of resolved point sources 
(such as active galactic nuclei (AGNs), pulsars, and star- 
forming galaxies), transient gamma-ray sources (e.g., 
gamma-ray bursts) , and the diffuse gamma-ray radiation 
comprised of emission from the Galaxy and the extra- 
galactic gamma-ray background (EGB). The origins of 
the EGB are, as yet, unknown; however, it is expected 
that unresolved, extragalactic point sources provide a siz- 
able contribution to the EGB. 

In its first year of taking data, Fermi observed 1451 
resolved point sources with significance g reater than 
5cr (F irst Fermi-LAT catalog (1FGL); lAbdo et al.l 
I2010H 1 of which 573 were associated with blazars 
(AGNs for which the jet is closely aligned with t he 
observer's line of sight (fBlandfor d fc Konigll [1979)). 
Thus, blazars constitute the largest class of astro- 
physical objects associated with gamma-ray sources, 
and it has long been suspected that unresolved 
blazars should provide a sub s tantia l contribution 
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It should then, perhaps, come as no surprise that 
the Fermi-LAT , with its improved sensitivity enabling 
it to observe many more blazars, measured the inte- 
grated intensity of the EGB to be substantially lower 
than that measured by the Energetic Gamma-Ray Ex- 
periment Telescope (EGRET) aboard the Compton 
Gamma-ray Observatory in the 1990s (1.03 (± 0.17) x 
10~ 5 cm -2 s _1 sr" 1 for Fermi versus 1.45 (± 0.05) x 
10~ 5 cm" 2 s- 1 sr" 1 for EGRET: ISreekumar et aDll99l 



lAbdo et"aDl2010gh . However, exactly how much of the 
change in the EGB integrated intensity is due to the 
ability to resolve many mo re blazars remains unc lear. 
On the one hand, since the ISreekumar et all ([1998D de- 
termination of the EGRET EGB, models of the galac- 
tic foreground emission have been updated to reflect re- 
cent o bservations of the interstellar medium (|Abdo et all 
l2010gf ) . accounting for at least some of the change. On 
the other hand, the determination of the distribution 
of blazars with respect to luminosity and redshift (the 
blazar gamma-ray luminosity function, GLF) and by ex- 
tension, their contribution to the EGB directly from their 
observed flux dist ribution (such as that performed in 
lAbdo et alll2lH0ch is non-trivial, for four reasons. 

First, blazars are variable at gamma-ray energies 
and without knowledge of the degree to which the 
gamma-ray luminosity of a blazar changes during 
the flaring state and the blazar duty cycle, the 
reconciliation of the observed bl azar source counts 
with model GLFs is impossible (Stcckcr fc Salamonl 
119961 ). This is because the increase in flux of a 
blaza r (sometimes up to an o rder of magnitude; see 
e.g.. iMcLaughlin. et all 119961: iMukheriee et all 119971: 
iNolan et al.ll2003t iVercellone et alll2004ft during the flar- 



2 



Venters & Pavlidou 



ing period introduces a selection effect in that flaring 
blazars are easier to detect than quiescent blazars, and 
quiescent blazars would have to be relatively bright and 
nearby to be detectable. The largest impact on the ob- 
served blazar source counts will be at fluxes just above 
the Fermi-hAT sensitivity as some blazars below the 
threshold make it into the sample because they are flar- 
ing. While the magnitude of the impact on the observed 
blazar source counts depends on the blazar GLF and 
the blazar duty cycle, the effect will be to flatten the 
faint-end slope of the observed blazar source counts (as 
blazars move from lower flux bins to higher ones) lead- 
ing to an underestimation of the blazar contribution to 
the E GB as extrapolated f rom the observed source counts 
alone. lAbdo et"aLl ()2010cD concluded that since the peak- 
to-mean flux ratio for most Fermi sources is ~ a factor 
of two, that there are no large systematic uncertainties 
due to blazar variability. However, a small peak-to-mean 
ratio is expected for blazars that are observed mostly in 
a flaring state, and as most of the observed blazars are 
fainter blazars that will be more represented by flaring 
blazars, then while it will be true that most blazars will 
have a small peak-to-mean ratio, it is not necessarily the 
case that the effect of blazar variability will be small. 

Second, the large angular resolution of the Fermi-hAT 
at lower energies (~ 5° at 100 MeV) could introduce sig- 
nificant source confusion for GLF models that predict 
high blazar densities resulting in fewer resolved blazars 
and a higher cont r ibution to the EGB at lower ener- 
gies ([Steckerl [l999t iStecker fc Venters! l20lTh . The im- 
pact of source confusion on observed source counts de- 
pends on the source density as given by the blazar GLF. 
As such, accounting for source confusion requires a pri- 
ori knowledge of the source density, which is exactly 
what the observer wishes to determine. Thus, source 
confusion further complicates analyses of the blazar con- 
tribution to the EGB based solely on observed source 
counts (for a detailed discussion of source confusion and 
the difference between the defin ition employed in this 
discussion and that emp loyed by lAbdo et aLil2010cl sec 
IStecker fc Ve^t^l2lLUl ). 

Third, as demonstrated in lAbdo et al.l (|2010cP ). the re- 
constructed fluxes of the individual sources are subject to 
a considerable amount of uncertainty, particularly at the 
faint end (see their Figure 6). Finally, the model of the 
Fermi- LAT detection efficiency employed by lAbdo et al.l 
(|2010cf ) assumed power-law spectra for the source WAs 
such, these uncertainties complicate the determination 
of the distribution of blazars with respect to luminos- 
ity and redshift (the GLF) and the determination of the 
blazar contribution to t he EGB from their observ ed flux 
distribution. Notably, IStecker fc Ventera (I2011P) foun d 
that contrary to the conclusion of lAbdo et al.1 (|2010d ). 
the observed flux distribution of Fermi blazars does not, 
as yet, rule out the possibility that the EGB is dom- 
inated by emission fr om unresolved blazars (see also 
lAbazaiian et al.l 120111) . Thus, the contribution of still 
unresolved blazars to the EGB remains in dispute. 

Additional information about the contributions to 

5 It should also be noted that the model for Fermi-hAT detec- 
tion efficiency (particularly at energies below ~ 300 MeV where 
multiple scattering in the detector becomes important) is subject 
to change during the course of observations by Fermi-LAT. 



the EGB can be obtained through studying the 
shape of its energy spectrum, and the shapes of the 
collective intensity spect ra of suspected contributors 
(IPavlidou fc Ventersl20"0l hereafter PV08 : 1 Venters et al.1 
2009; Vc ntersll20101 ). Analysis of the first year of Fermi 
data yielded an EGB spectrum consistent wi th a feature- 
less po wer law with spectral index V ~ 2.4 ()Abdo et al.l 
2010g). Upon first reflection, the Fermi EGB spectrum 
appears to be consistent with the hypothesis that the 
EGB is dominated by emission from unresolved blazars 
since the mean spectral index for blazars, To, is also 
~ 2.4. However, several effects complicate this simple 
picture. 

First, the astrophysical population of blazars is ac- 
tually composed of two separate sub-populations (flat- 
spectrum radio quasars, FSRQs, and BL Lacs) with 
distinct spectral prop erties (r n ~ 2.45 f or FSRQs and 
Tq ~ 2.2 for BL Lacs; lAbdo et al.ll2010d ). Second, even 
within a given sub-population of blazars, the spectral in- 
dices of i ndividual blazars fo r m a d i stribution with some 
sprea d (jStecker fc Salamonl 119961 : iVenters fc Pavlidou] 
l2007h . which causes the collective spectrum of unre- 
solved blazars to curve as harder bl azar spectra become 
more important at higher energ ies (Stcckcr fc Salamonl 
119961: IPavlidou fc Venters! [20081. On the other hand, 
this effect is somewhat mitigated by the spectral bias 
introduced by the fact that blazars with harder spec- 
tral indices are easier to observe in a flux-limited sur- 
vej@, and thus, are more likely to be resolvable and not 
play as big a role in producing the EGB. Third, just 
as unresolved blazars are expected to contribute to the 
EGB, unresolved members of other known astrophysi- 
cal gamma-ray emitters (such as star-forming and star- 
burst galaxies) should also contribute, but with spec- 
tra that are substa ntially different from those of blazars 
(IPavlidou fc Fields! ]2002t iFields et "all 120101: [Lacki et al.l 
120111: iMakiva et al.l 120111: IStecker fc Venters! 1201 ID and 
may not even resemble power laws at gamma-ray ener- 
gies. Finally, recent observations conducted by Fermi in- 
dicate that even blazar spectra can break at ~ GeV ener- 
gies a nd would no longer be describable by s imple power 
laws (lAbdo et al.1 l2009ab llnoue et al.l 120081: llnoue et al.1 
[20101: llnoue fc Tot anil 1200^ 

The intuitive conclusion would be that the combina- 
tion of these effects should cause the energy spectrum of 
the E GB intensity to exhibit features (jPavlidou fc Fields! 
2002) . The observation instead of an EGB with a single, 
featureless power-law spectrum begs the question: is the 
lack of observed spectral features necessarily indicative of 
a single dominant source population with individual un- 
broken power-law spectra in the energy range of > 300 
MeV? If not, then what can such a featureless power- 
law EGB spectrum tell us about the relative contributions 
of the separate populations and how do they depend on 
their collective spectra? The answer to these questions 
provide novel constraints on blazar GLFs and thus could 
have profound implications for the cosmological prop- 

6 It should be noted that the Fermi survey is not exactly a 
flux-limited survey due to the non-uniformity of the total diffuse 
background throughout the sky. As such, spectral bias is more sig- 
nificant in the Fermi survey than typical for a truly flux-limited 
survey. However, as evid enced from the anal ysis of the mean spec- 
tral indices performed by Abdo et al. (2010c), this appears to have 
more of an effect on BL Lacs than FSRQs. 
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erties of known astrophysical gamma-ray emitters, their 
corresponding contributions to the EGB, and the general 
properties of their spectra at gamma-ray energies. 

In this paper, we calculate the collective spectrum 
of unresolved blazars with individual spectra exhibiting 
broken power laws. In so doing, we seek to determine 
whether such a population of unresolved blazars can re- 
sult in a featureless power-law spectrum resembling that 
of the EGB. If so, we could then investigate the impli- 
cations for the relative contributions of FSRQs and BL 
Lacs, which could, in turn, have implications for the cos- 
mological properties of blazars. In Section 2, we present 
the formalism of the calculation of the spectral shape of 
the collective unresolved blazar emission. In Section 3, 
we discuss the inputs of the calculation and their un- 
certainties. In Section 4, we present the results of the 
calculation, and we discuss them in Section 5. 

2. FORMALISM 

To calculate the collective spectrum of unre solved 
blazars, we follow the formalism as outlined in IPV08I 
with one major difference. Instead of taking blazar 
spectra to be simple power laws over gamma-ray ener- 
gies (F oc E~ r , where T is the photon spectral index 
at gamma-ray energies), blazar spectra are taken to be 
smoothly broken power laws: 



Fe(E) — Fq 



E 
E~ b 



rin 



r 2 n 



-1/n 



(1) 



where Fe(E) is the differential photon flux in units of 
photons per unit area per unit energy per unit time, Eb 
is the break energy, is the low-energy slope, T 2 is 
the high-energy slope, and n quantifies the sharpness of 
the transition from the low-energy power law to the high- 
energy power law. For the purposes of this paper, we take 
n=l. The total flux, F, of photons with energies greater 
than some fiducial energy, Ef, is found by integrating 
F E (E) over energy, 



fin 



E 
E, 
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dE. 



(2) 



For the purposes of this paper, we take Ej — 100 MeV. 
Then, the contribution of a single unresolved blazar to 
the EGB is 



h 



F [{E/E b ) T ^ + [E/E b ) T - n ] 1/n 
4tt JT [{E/E b y^ + (E/E b )^]- 1/n dE 



(3) 



where the flux of one source is uniformly distributed 
over the entire sky in anticipation of an isotropically dis- 
tributed cosmological population and I has units of pho- 
tons per unit area per unit energy per unit time per unit 

solid angle. 

Following IPVOl we characterize the flux distribution 
of unresolved blazars as a function g(F) and the distri- 
bution of blazar spectral indices (or spectral index dis- 
tribution, SID) as a function p(T), where Ti = T — AFi, 
T 2 = T + Ar 2 , and p(r), Ari, and Ar 2 are determined 
from observations (see Section [3]). Then, the total con- 



tribution of unresolved blazars to the EGB is given by 
dT / dFg{F)I lP {T) , (4) 

-oc JO 

where -Fmin is the sensitivity of the gamma-ray telescope 
under consideration. For the first year of Fermi data, 
-Fmin ~2x 10~ 9 photons cm -2 s -1 . 

Equation (U) can be characterized in terms of factors 
that determine the overall magnitude of the unresolved 
blazar contribution to the EGB (flux terms) and factors 
that determine the overall shape of the blazar contribu- 
tion (spectral index terms). For a carefully chosen defi- 
nition of p(T) (determined from the analysis of spectral 
indices of the flux-limited sample of Fermi blazars that 
accounts for the spectral bias inherent in a fl ux-limited 
catalog; see Section [3l and [Venters et aLll2009f) . the mag- 
nitude and shape terms decouple, and Equation (j4|) can 
be rewritten as 



1(E) = 

p OC 

Jo dTp{T) 



[(E/E b )( r - Ar ^ n + (E/E b )( r + Ar ^ 
S(E f ,T) 



-l/r, 
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where Iq is a normalization constant depending on the 
flux distribution of unresolved blazars and 
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3. INPUTS 

As demonstrated in IPV081 the unresolved blazar con- 
tribution to the EGB is not just a question of magni- 
tude, but also of spectral shape, and the spectral shape 
is sensitive to the distri bution of bl azar spectra l indices 
at GeV energies. Both |PV08l and lAbdo et all (|2010d ) 
assumed that blazar spectra at gamma-ray energies are 
power laws and did not account for possible breaks in the 
spectra. In this paper, blazar spectra take the forms of 
smoothly-broken power laws (see Section [2J . 

For each sub-population of blazars, we deter- 
mi ne the SID, jp(r), from the likelihood analysis 
of iVenters fc Pavlidoul (|2007l ) fittin g blazars from the 
Fermi-LAT First Catalog of AGNs (|Abdo et al.ll2010dl ) 
to Gaussian SIDs accounting for errors in measurement 
of individual spectral indices. Due to the survey bias 
towards harder spectral indices present in Fermi data, 
we applied the likelihood analysis only to the subset of 
blazars with photon fluxes > 7 x 10~ 8 photons cm" 2 s" 1 
and galactic latitudes > 10° (as per lAbdo et alJl20l0ch . 
We determined that the maximum-likelihood Gaussian 
SID can be characterized by a mean (ro) and a spread 
(do) with maximum- likelihood parameters determined to 
be r = 2.45 and <t = 0.16 for FSRQs and T„ = 2.17 
and do = 0.23 for BL Lacs, which are similar to the find- 
ings of lAbdo et"ail (|2010dl . Based on 1FGL spectra of 
two pro minent blazars from the Fermi-LAT Bright AGN 
Sample (|Abdo et al.ll2009al ; see Figure [T]), we model the 
spectral breaks by taking AFi =0.1 and Ar 2 = 0.9. As 
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in the measured blazar spectra, we treat the break en- 
ergies distinctly for FSRQs and BL Lacs. We took the 
break energy, E b , to be 4 GeV for FSRQs and 15 GeV 
for BL Lacs. In so doing, we consider two cases: 

Scenario 1. - Blazars within a population evolve 
such that their break energies are observed to be 
roughly the same. 

Scenario 2. - The break energies of the intrinsic 
spectra of blazars within a subpopulation are the 
same, but because of redshift effects, they are ob- 
served to be different. 

In order to perform the calculation in the Scenario 2, 
we separately model the contribution arising from dif- 
ferent redshift bins from an assumed GLF model for a 
given subpopulation of blazars and then determine the 
composite spectrum in the final step. For the purposes of 
this analysis, we assume the best-fit p u re lum inosity evo- 
lution model of iNarumoto fc Totanil (|2006t l for FSRQs 
and their best-fit luminosity-dependent density evolution 
model for BL Lacfl In this manner, we seek to inves- 
tigate the impact of changing the GLF model(s) of the 
blazars. We should also note that alternati ve statistical 
treatm ents on more extensive data sets as in lAbdo et al.1 
(|2010ff ) could result in different estimates of the break 
energies and changes in spectral index. 

Another alterna tive treatment performed by 
lAbdo et"ail (|2010d) consists of stacking the spectra 
of observed blazars. However, the analysis was per- 
formed on the flux-limited sample of Fermi blazars, 
which consists almost entirely of FSRQs (hence, the sim- 
ilarity between the stacked spectrum for the flux-limited 
sample of blazars and that of FSRQs). Moreover, while 
such an analysis is effective in determining the average 
spectral index of a given population of blazars, it reveals 
nothing about the spectral properties of the population 
of blazars (e.g., the spread in the SID, spectral breaks). 
Furthermore, as the Fermi-LAT survey is biased against 
sources with softer spectra, it is certainly biased against 
sources exhibiting significant breaks in their spectra. 
We also note that despite findin g that the S IDs o f 
blazars have non-zero spreads, lAbdo et all (|2010d l 
did not account for these spreads in their calculated 
spectra of the blazar contributions to the EGB as 
evidenced by the fact that these spectra, as well as 
their estimated uncertai nties, are still powe r laws (see 
Figures 18, 19, and 20 of lAbdo et aTll2010d ). Similarly, 
lAbdo et"aT1 (j2010d ) did not account for spectral breaks 
in their calculated spectra of the blazar contributions 
to the EGB. Inclusion of these effects would introduce 
curvature in the spectra of the blazar contributions to 
the EGB and not r esult in the power laws indicated by 
lAbdo et all (l20T0cT) . 

As in PV08, we do not include information about the 
magnitude of the unresolved blazar contribution to the 
EGB. Instead, we normalize the collective spectrum of 
each sub-population in order to best fit (as determined 
using a x 2 -analysis) the measured s pectrum of the EGB 
from the Fermi-LAT first year data (|Abdo et al.ll2010gf ). 
For the composite spectrum of FSRQs and BL Lacs, we 

7 In so doing, we account for the observation that BL Lacs are 
on average situated at lower redshifts than FSRQs. 
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Fig. 1. — Sample broken power-law spectra for blazars. Top: 
PKS0235+164 measured spectrum with broken power-law (solid) 
and single power-law spectra (dashed; T = 2.1433). Bottom: 
3C454.3 measured spectrum with broken power-law and single 
power-law spectra (r = 2.4662). 

normalize the collective spectrum of each sub-population 
such that their total composite spectrum is the best-fit 
spectrum to the measured spectrum. 

We should note that while applying the aforemen- 
tioned cuts to t he sample accounts for much (though 
likely not all; see I Venters et al.l l2009') of the survey spec- 
tral bias, doing so also significantly reduces the sample 
size and will ultimately introduce more uncertainty in 
the determination of the likelihood parameters. Further- 
more, while in the flux-limited sample of blazars, spec- 
tr al index and flux d o appear to decouple (see Figure 1 
of lAbdo et al.l I20"l0cf) . we expect it to be only approx- 
imately correct. A more detailed calculation, however, 
requires a statistical analysis to determine the intrinsic 
blazar SID, which would, in turn, require knowledge of 
the blazar GLFs since the survey spectral bias depends 
on the redshift and luminosity distributions of blazart0 

8 For this same reason, we also neglect source confusion at low 
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(| Venters et al.l l2009f ). In the case of FSR Qs, applying 
the high-flux cut changes the S ID very little (jAbdo et al.l 
l20lMIStecker fc VentersllMll) . so we do not expect the 
collective spectrum of unresolved FSRQs to be appre- 
ciably different from that determined from the resolved 
FSRQs. On the other hand, in the case of BL Lacs, the 
sparseness of the population at fluxes above the cut in- 
troduces a considerable margin of error. As such, the col- 
lective spectrum of unresolved BL Lacs could be harder 
or softer than that expected from the resolved BL Lacs. 
As BL Lacs are, on average, harder than FSRQs, they 
will likely contribute most significantly at the higher en- 
ergies, and hence, the differences in spectra between re- 
solved and unresolved BL Lacs will mostly impact the 
collective unresolved blazar spectrum at higher energies. 

Finally, we should note that recent multi-wavelength 
observations conducted by Fermi and other telescopes 
suggest that blazars can be further subdivided into the 
categories of low synchrotron peak, inter mediate syn- 
chrotr on peak, and high synchrotron peak ()Abdo et al.l 
l2010el ). In principle, one could gain further insight by 
applying this procedure to each subdivision as there is 
some indication that these cate gories are also spec trally 
distinct at gamma-ray en ergies (lAbdo et alj l2010en. No- 
tably, as demonstrated in lAbdo et al.l (|20i*oi ) . most high- 
synchrotron peaked BL Lacs (HBLs) do not exhibit spec- 
tral breaks in Fermi-LAT energy range while those HBL 
breaks that have been observed are different in character 
(going from soft to hard rather than hard to soft). How- 
ever, in light of the cuts on galactic latitude and flux, 
further subdividing the sample of blazars could lead to 
small sample sizes in some of the categories, and conse- 
quently to poor SID parameter determination. Thus, for 
the purposes of this analysis, we retain the original clas- 
sification scheme of FSRQs and BL Lacs. In any case, 
given that BL Lacs likely contribute most significantly 
at higher energies, the effect of the HBL spectra is likely 
to improve the fit at high energies, though a more de- 
tailed study is in order as more data become available. 
One might also be concerned about systematic changes 
in blazar spectral indices while flaring; however, analy- 
ses of EGRET blazar spectral indices found no evidence 
of s uch systematic changes in spectral index with flar - 
ing (|Nandikotkur et al.ll2057t iVenters fc Pavlidoull2007l ). 
and Fermi observations of individual blazars have thus 
far revealed no systematic changes in spectral index with 
time or flux (|Abdo et al.ll2009Ul2010atlAckermann et al.l 

[SlGJ. 

4. RESULTS 

The spectral shapes of the contributions of unresolved 
blazars to the EGB as determined from the Gaussian 
SIDs and broken power-law gamma-ray spectra discussed 
in Section [3] are plotted in Figures HHU In Figure [5J we 
have plotted the best-fit spectral shapes of the collec- 
tive intensities of unresolved FSRQs (x^cd = 4-1) and 
BL Lacs (Xred = 3-3). In Figure El we have plotted the 
best-fit spectral shape of the total collective spectrum of 
unresolved blazars in Scenario I (solid line; Xrcd = 0.7; 
for explanation, see Section [3]) and the individual con- 
tributions from FSRQs (dashed line) and BL Lacs (dot- 
energies arising from the increase in Fermi-hAT point-spread func- 
tion below 1 GeV. 
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Fig. 2. — Best-fit shapes of the individual collective intensities of 
the sub-populations of unresolved blazars. Top: the best-fit shape 
for FSRQs. Bottom: the best-fit shape for BL Lacs. Data points: 
spectr um of the EGB as measured by the Fermi-hAT f Al b do et al.l 
2010a). 



dashed line). In Figure HI we have plotted the best- fit 
spectral shape of the total collective spectrum of unre- 
solved blazars in Scenario 2 (x?od = 0.6). For compar- 
ison, we have also plotted the best-fit spectral shape of 
the total collective spectrum assuming that FSRQ spec- 
tra break at Eb ~ 1 GeV and BL Lac spectra break at 
Eb ~ 5 GeV (assuming Scenario 1; Figure [5]). Note that 
such spectra do not adequately fit the spectra in Figure 
[TJ As such, this alternative is included merely for the 
sake of comparison. 

As noted in Section [3l BL Lacs tend to be harder than 
FSRQs, and thus, the BL Lac collective spectrum is no- 
ticeably harder than that of the FSRQs. Also apparent 
is the effect of the higher break energy for the BL Lacs 
as compared with the FSRQs. As can be seen in Figure 
[H the collective spectrum of neither sub-population of 
blazars reproduces well the spectrum of the EGB: the BL 
Lacs are too hard while the FSRQs are too soft. How- 
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Fig. 3. — Best-fit shape of the total collective intensity of unre- 
solved blazars in Scenario 1 as described in Section[3] Dashed line: 
the contribution from FSRQs. Dot-dashed line: the contribution 
from BL Lacs. Solid line: the shape of the combined population of 
FSRQs and BL Lacs. Data points: same as in Figure [2] 




10 



10 10 
Energy [GeV] 



10 



Fig. 4. — Same as in Figure[3] but for Scenario 2 as described in 
Section \3\ 



ever, if the spectra are added together as in Figures |3J 
and 0] (renormalized as described in Section |3j, the re- 
sulting spectrum (solid line) reproduces well that of the 
EGB (being within the statistical error bars of nearly all 
of the data points and within ~ 1.5a of the last data 
point). The harder spectra and higher break energies for 
the BL Lacs compensate for the suppressed intensity of 
FSRQs at higher energies. Notably, the redshift depen- 
dence of Scenario 2 shifts the breaks to lower energies 
resulting in the transition from FSRQ dominance to BL 
Lac dominance occurring at lower energies (~ 1 GeV 
rather than ~ 3 GeV). Thus, the relative contribution 
from BL Lacs to the total blazar collective spectrum is 




10 10 10 10 

Energy [GeV] 

Fig. 5. — Same as in Figure [3] but for lower break energies. 



greater in Scenario 2 than in Scenario 1. In the case of the 
even lower break energies as shown in Figure [5j the total 
collective spectrum fits the EGB fairly well (x;?ed = 1-1) j 
but the relative contributions of the FSRQs and BL Lacs 
are such that the contribution from BL Lacs dominates 
at all energies. This is in contrast with the scenarios 
presented in Figures |3] and U in which BL Lacs dominate 
only energies greater than ~ few GeV. Also, the sharper 
contrast between Figures [3] and \5\ than that between Fig- 
ures [3] and [5] indicate that the transition and the resulting 
relative contributions of FSRQs and BL Lacs are more 
sensitive to the break energies than model of the blazar 
GLF(s). This is due to the fact that each contribution 
is dominated by its closest and brightest (though still 
unresolved) members for which the observed break ener- 
gies are similar to their intrinsic values for the GLF(s) 
considered. 

5. DISCUSSION AND CONCLUSION 

We have calculated the spectral shapes of the contri- 
butions to the EGB arising from BL Lacs and FSRQs as- 
suming blazar spectra can be described as broken power 
laws. We found that in the case that blazar spectral 
breaks are indeed common neither sub-population alone 
can adequately reproduce the spectrum of the EGB. 
However, in a combined spectrum, the harder spectra of 
the BL Lacs could compensate for the softer spectra of 
the FSRQs, resulting in a collective blazar spectrum that 
is similar to that of the EGB within uncertainties. Fur- 
thermore, we have found that the relative contributions 
of FSRQs and BL Lacs required to fit the EGB spectrum 
are sensitive to the nature of the spectral breaks and, to 
a lesser extent, the blazar GLFs. 

The question of whether the relative contributions nec- 
essary to reproduce the overall spectrum of the EGB as 
determined in this method are reasonable given future 
Fermi measurements of the blazar GLFs (which account 
for the notable observational uncertainties discussed be- 
low and in Section Q]) would provide insight into the over- 
all blazar contribution to the EGB as well as the nature 
of blazar variability. According to current Fermi obser- 
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vations, the luminosities of FSRQs ar e greater than those 
of BL Lacs (L 7 ,fsrq ~ 10 2 x £ 7 ,bll; lAbdo et ai1l20Q9al 
2010d). On the other hand, BL Lacs are li kely situated a t 
lower, on average, redshifts than FSRQs (|Dermerll2007h . 
Thus, it is possible that the closer proximity of BL Lacs 
(or their numbers) compensates for their deficit in lumi- 
nosity with respect to FSRQs resulting in their roughly 
comparable or slightly enhanced relative contribution to 
the collective blazar spectrum seen in Figures 02 [4] and 
O However, it should be noted that the closer proximity 
of BL Lacs and their harder spectral indices would also 
make them more easily observable by Fermi, limiting the 
contribution of unresolved BL Lacs to the EGB-though 
if the luminosity function for BL Lacs is broad enough, 
there could be many low-luminosity BL Lacs that would 
escape detection. 

It thus remains to be seen how the collective intensi- 
ties of the sub-populations of blazars actually compare 
with one another. In order for a true comparison to be 
drawn, the GLF for each sub-population of blazars needs 
to be measured. If it should be the case that the relative 
contributions are not reasonable given gamma-ray ob- 
servations, then either the breaks included in the best-fit 
spectrum are not typical of blazars or blazars are not suf- 
ficient to explain the EGB. Already, Fermi studies of ob- 
served blazar flux distributions have suggested that emis- 
sion from unresolved blazars m ay not comprise th e domi- 
nant contribution to the EGB (jAbdo et al.ll2010cl ). How- 
ever, given that the physics behind the gamma-ray emis- 
sion of other known and speculated contributors (i.e., 
normal galaxies, cascades of ultra-high energy cosmic 
rays and TeV gamma rays, and dark matter annihila- 
tion) result in spectra that are quite distinc t from that 
of the EGB as measur ed by Fermi (see e.g.. | Ando et aLl 
2007t iKalashev et all l2009t iSiegal-Gaskins fc Pavlidoul 
2009tlAhlers et al.ll2010|: iFields et al.ll2010tlVentersll201d : 
Berezinskv et al.ll2011t iStecker fc Venters! [20 111) , a close 
resemblance of the collective spectrum of unresolved 
blazars could be striking. 

The necessary reconciliation of the clues to the blazar 
contribution to the EGB provided by studies of the col- 
lective blazar spectrum with the clues provided by stud- 
ies of the blazar source counts could thus provide insight 
in the gamma-ray emission properties of blazars. Blazar 
variability plays a substantial role as flaring blazars 
would be more easily observed by Fermi than quiescent 
blazars. As such, studies of the blazar GLF, observed 
blazar flux distributions, and ultimately, the blazar con- 
tribution to the EGB are largely a question of the blazar 
duty cycle. If blazars spend the majority of their lifetimes 
in the more observationally challenging quiescent state, 
then studies of the blazar contribution to the EGB based 
on observed flux di stributions such as that presented in 
lAbdo et al.l ()2010cf ) could underestimate the number of 
quiescent blazars resulting in an underestimation of the 
blazar contribution to the EGB. Thus, the apparent dis- 
crepancy of the predictions of the two analyses could be 
the result of the variability of gamma-ray emission in 
blazars, the study of which could have implications for 
gamma-ray emission in blazars. 

The comparison of the two types of analyses could also 
provide insight in the cosmological properties of blazars. 
At 100 MeV, the Fe rmi-LAT angular resolution is ~ 5° 
(jAtwood et al.l [20091) ; hence, for blazar GLFs that pre- 



dict large number densities, many blazars, particularly 
those on the faint-end of the source count distributions, 
that are, in principle, observable would not be distin- 
guishable from other blazars (and thus are unresolved) 
and would not be included in source count distributions. 
In such a scenario, the faint-end of the blazar flux distri- 
bution would underestimate the number of blazars in a 
given flux bin, and the faint-end slope might appear flat- 
ter than it should bc0, resulting in an underestimation of 
the blazar contribution to the EGB. Thus, the apparent 
discrepancy of the results of the two analyses could be 
the result of a large number density of blazars. 

In conclusion, we have demonstrated that even with 
the inclusion of spectral breaks, the collective spectrum 
of unresolved blazars reproduces well the spectrum of the 
Fermi EGB for several break models. As such, we find 
that the possibility that the collective intensity of un- 
resolved blazars dominates the EGB is not excluded on 
spectral grounds, even if spectral breaks are shown to be 
common among Fermi blazars. Given the remaining con- 
troversy concerning the blazar GLFs, we conclude that 
it is, as yet, premature to rule out blazar dominance of 
the EGB. Furthermore, we have shown that relative con- 
tributions of the sub-populations of blazars required to 
fit the EGB spectrum are sensitive to the nature of their 
breaks; hence, the methodology we present in this paper 
can be used to constrain the GLFs of blazars. As models 
for the GLFs of blazars (accounting for the uncertain- 
ties outlined in this paper) and more data on spectral 
breaks become available, the study of the spectral shape 
of the blazar contribution to the EGB in light of the 
breaks can provide insight into the high-energy physics 
of blazars and their cosmological properties. 
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9 Another symptom (but not proof) of the effect of source con- 
fusion would be a similar Fermi measurement of the EGB as that 
of EGRET at lower energies since at these energies, the angu- 
lar resolution of Fermi - LAT is compa r able t o that of EGRET 
IIThompson et alj [19931 ; lAtwood etaTl [20091: IStecker fc Venters! 
120111) . Intriguingly, the first few d ata points of t he Fe rmi EGB 
do appear to be consistent with the Strong et al. (2004) determi- 
nation of the EGRET EGB. 
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